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Abstract. We consider isotropic version of XY model in transverse magnetic field 
[in one space dimension]. We calculated of asymptotic of temperature correlations. 
We represent the quantum correlation function as a tau function of a classical 
completely integrable differential equation. This is the well-known Ablowitz-Ladik 
lattice equation. 
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1. Introduction 

The XY model was introduced and studied by E. Lieb, T. Schultz and D. Mattis [1]. 
It describes the interaction of spins 1/2 situated on a 1-dimensional periodic lattice. 
The Hamiltonian of the model is 



n 

Here a are Pauli matrices, h is transverse magnetic field and n enumerates the sites of 
the lattice. At zero temperature the problem of evaluation of asymptotics of correlation 
functions was solved in [2,3]. Here we consider the temperature correlation function. 



g{n,t) ^ ^ , n = 712-711, ^ = ^2-^1 (2) 

Tre T 

for the infinite lattice. We consider finite temperature < T < 00 and a moderate 
magnetic field < /i < 2. 

We evaluated the asymptotics in cases where both space and time separation go 
to infinity 7i — > 00, i ^ 00, in some direction </? 

Ti TT 

— = cot(/?, 0<(^<-. (3) 

In accordance with our calculations, correlation function g{n, t) decays exponentially 
in any direction, but the rate of decay depends on the direction. In the space like 
direction, < 9? < ^, the asymptotics are 



/ \ ^ f ^ /"'^ 7 1 1 //i — 2cosp, 1 
g{n,t) ^Cex^l— J dp In tanh( ) j 

In the time like direction f < </? < f , the asymptotics are different: 



(4) 
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g{n, t) Ct^-^-l+-'--) expi^^J dp \n - At sin p\ In 



tanh( 



h — 2 cosp 
T 



(5) 



The values which define the pre-exponent, are 



— — In 

+ 27T 



tanh 



h — 2 cos po 
f 



i/_ = — In 
27r 



tanh 



h + 2 cos Po 
f 



(6) 



where ^ = sinpo- Equation (5) is vahd in the whole time like cone, with exception 
of one direction h = 2cospo- Higher asymptotic corrections will modify formulae by a 
factor of {l + c{t, x)) (c decays exponentially in the space-like region and as t~^^'^ in the 
time-like region). Also, it should be mentioned that the constant factor C in (4) does 
not depend on the direction ip, but does depend on ip in (5). We want to emphasize that 
for the pure time direction, (p — 7r/2, the leading factor in the asymptotics (exponent 
in (5)) was first obtained in [10]. 

To derive these formulae we went through a few steps. 

The first step: The explicit expression for eigenfunctions of the Hamiltonian (1) 
(see [1]) was used to represent the correlation function as a determinant of an integral 
operator (of Fredholm type) [8]. In order to explain we need to introduce some notation. 
Let us consider the integral operator V. Its kernel is equal to 



e+(A)e_(/i) - e_(A)e+(//) 



7r(A - //) 

Here A and /i are complex variables, which go along the circle |A| 



(8) 

1 in the 
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positive direction. The functions e± are 

e_ (A) = A-"/' • e-'*(^+^/^) (9) 



where 



^;(A) = <^ 1 + exp ^ ^' 



T 



(10) 



and 



e+(A) = e_(A)^(x,t,A). (11) 
Here E is defined as an integral 

E{n,t,\)^-v.p. [ exp{2ii(/x + -)}-^^. (12) 
n J IJ, IJ,- A 

It is convenient to define functions /±(A) as solutions of the following integral equations: 

iI + V)fk^ek. (13) 
Here / is the identity operator and k = ±. Next we define the potentials Bkj{k,j = ±): 

BkAn,t) = ^J MX)eAX)^. (M) 

They depend on space and time variables n, t. These we shall use to define new 
potentials b^j'. 

b—{n,t) = B__{n,t), 

(15) 

6++(n,t) = B++{n,t) - 2iG{n,t)B+_{n,t) - G{n,t). 



Temperature Correlation of Quantum Spins 



5 



Here we used the function 



G{n, t) = ^ [ A"-i exp{2it{X + \)}dX. (16) 
zm J A 

Now all the notation is ready to write a determinant formula for the correlation function 

g{n,t) (see (4)): 



g{n,t) = e-^"''b++{n,t)exp{a{n,t)}. (17) 
Here e'^ is a determinant of the integral operator 

exp{(7(n, t)} = det(l + V). (18) 

Second step: Formulae (8)-(15) can be used to show that the potentials 6++ and 
b satisfy a system of nonlinear differential equations. 

^-6__(n, t) = (1 + 46__(n, t)b++{n, t)){b..{n + l,t) + b..{n - 1, t)) 



-'--b+^{n, t)^{l + 46__(n, t)b++{n, t)){b++{n +l,t)+ b++{n - 1, t)). (19) 

The derivation of these equations is similar to [4,5,7]. Equations (19) arc completely 
integrable differential equations. They were first discovered by Ablowitz and Ladik [9] 
as an integrable discretization of the nonlinear Schroedinger equation. The logarithmic 
derivatives of a{x,t) (see (18)) can be expressed in terms of solutions of the system 
(19): 
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= 26__(n, t)b++{n, t) - b++{n - 1, t)h^^{n + 1, t)- 

-6__(n- l,t)6++(n+ - 46++(n, t)6__(n, t)[6++(n - 1, t)6__(n + 1, t)+ 

l,i)] 

(20) 

cr(n + 1, t) + - 1, - 2(j(n, = ln[l + 46__(n, ^)] (21) 



— [(T(n + 1, i) - (7(n, i)] = 8i[6++(n + 1, i)&__(n, i) - 6++(n, i)6__(n + 1, t)]. (22) 

This shows that the quantum correlation function g (2) can be expressed in terms 
of the solution of the system (19). The meaning of all these formulae is that the 
correlation function of the XY model is the r- function (in a sense of the well-known 
works [11,12]) of Ablowitz-Ladik's differential-difference equations. In the papers [5,6,7] 
the relation between the r functions of the classical partial differential equations and 
quantum correlation functions, together with the history of the question, is explained in 
more detail. It is also worth mentioning that the idea to connect quantum correlation 
functions and classical completely integrable systems goes back to the work [13] and 
was first apphed to the XY model in [14]. 

Third step: In order to evaluate the asymptotics one should solve Ablowitz-Ladik's 
differential equation. Initial data can be extracted from the integral representations 
(8)- (18). We use the Riemann-Hilbert problem in order to evaluate the asymptotics of 
the solution of equation (19). It is quite similar to the nonlinear Shrodinger case [6,7]. 

Finally, let us explain the physical meaning of our asymptotic formula (4). We 
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start from the expression for the free energy [1]: 

T r , . Acosp-2h,\ 
f{h) = -h-—Jjpln\^l + exp[ 1 ]j. (23) 

We emphasize the dependence on the magnetic field h. The definition of f{h) is 
standard: 

Tre-f = eM~fm- (24) 

Here L is the length of the box. Let us use Jordan- Wigner transformation to transform 
correlator (2) (in the equal time case) : 



n2-l ^ ^ 

<(0)<7-(0) = Vn.exp<; z^ ^l^t^k^-il-al). (25) 

fc=ni+l J 

Here i/jk is a canonical Fermi field. We note that numerator in (2) differs from the 
denominator by replacement of the magnetic field h ^ h — i7rT/2 on the space interval 
[ni + 1, ^2 — 1]. This leads us to the following asymptotic expression for correlator 

g{n,0) 



g{n, 0) exp < Re 



n 
'f 



f{h) - f{h - 



inT , 



(26) 



The reason we wrote Re is that ivr in (25) can be replaced by — ivr. It is remarkable 
that (26) coincides with the correct answer (4). It is also worth mentioning that to 
go to the exponent in (5) one should replace the differential d{np) by the expression 
\d{np — te{p)\, where e{p) — — 4cosp + 2h is the energy of the quasiparticle of the 
model. 
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